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ABSTRACT 

We propose a nonhnear a-model in a curved space as a general in- 
tegrable elhptic model. We construct its exact solutions and obtain 
energy estimates near the critical point. We consider the Pohlmeyer 
transformation in Euclidean space and investigate the gauge equiva- 
lence conditions for abroad class of elliptic equations. We develop the 
inverse scattering transform method for the sinh-Gordon equation and 
evaluate its exact and asymptotic solutions 
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1. INTODUCTION. 



This work is devoted to studying a broad class of problems involving two-dimensional 
nonlinear integrable elliptic equations. These equations, which describe stationary pro- 
cesses in nonlinear media, usually, govern the spatial distributions of physical fields 
whose value are fixed at the boundary of the domain under consideration. The re- 
sulting boundary- value problems have a characteristic feature: as the "observer" pen- 
etrates deeper into the medium (which means considering the asymptotic behavior as 
r"^ = + y"^ ^ oo), the field in question become weak, and the equations become the 
classical linear equations of mathematical physics. For this reason, the inverse scattering 
transform method (ISTM) applied to the elliptic equations (as well as the hyperbolic 
ones) is a nonlinear analogue of the Fourier transform method, a tool often applied to 
linear problems. This analogy allows one, in particular, to interpret the discrete-spectrum 
solutions of nonlinear elliptic equations as " solitonlike" excitations (vortices, or defects) 
in nonlinear condensed media. These excitations possess some properties of solitons and 
also some special properties of their own. 

In [1-4] , the boundary- value problems on the half-plane were posed for two-dimensional 
elliptic equations, the direct and the inverse scattering problems were investigated for the 
corresponding operators of the associated linear problems, and the exact solutions, the 
conservation laws, and the trace identities were also found. In addition, the boundary- 
value problems for the elliptic versions of the two-dimensional Heisenberg ferromagnet 
(of the 0(3) a- model) were shown to be gauge equivalent to the sinh-Gordon equation 
[4]- [5]. These results have thus allowed us to extend the methods used with the hyperbolic 
equations to the elliptic case, which has numerous important physical applications. 

At the same time, a currently important problem consists in constructing more-realistic 
models that would still be integrable on the one hand and would account for various 
features of the physical systems on the other hand. Solving one such problem was begun 
in [6], where a cr-model in a curved space was proposed (a magnet with a variable nominal 
magnetization). We show in what follows that this model (more precisely, its integrable 
version) is the most general of the presently known physical models that are described by 
elliptic equations and possess Lax representations. 

In this paper, which is an expanded version of [6, 1] and a revised version of [7], our 
aim is to obtain exact solutions of this model, to investigate its properties, and also to 
find a chain of transitions leading from this model to those that have been known and 
investigated before. For one such model (sinh-Gordon), we use the ISTM to construct 
exact solutions and give their physical interpretation. 

The paper is organized as follows. In Sec. 2, we consider a cx-model in a curved space, 
find its exact solutions, and investigate some thermodynamic properties. In Sec. 3, we 
propose a Euclidean version of the Pohlmeyer transformation and investigate the question 
of the gauge equivalence (GE) of a broad family of boundary-value problems. In Sec. 4, 
we develop the ISTM for a model of charged particles on a half-plane, find a series of its 
exact solutions, and also find the conservation laws and trace identities. 
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2. NONLINEAR a-MODEL IN CURVED SPACE. 



As is well known, the condition that the absolute value of the magnetic moment 
be constant is an essential assumption in the phenomenological description of weakly 
excited states of magnets. With this assumption, the evolution of system states amounts 
to the rotation of the magnetization vector. A similar statement is true for magnet 
models described by the Landau-Lifshits equation. Such models provide a physically 
adequate picture far from the Curie point. For moderate temperatures, however, one can 
expect a spatial variation of the absolute (nominal) magnetization of the sample material. 
This hypothesis allows one to investigate the magnet near the critical point (but outside 
the fluctuational region) using the phenomenological approach and to describe a smooth 
spatial transition to the paramagnetic phase. 

We consider the two-dimensional isotropic Heisenberg ferromagnet, i.e., the two- 
dimensional stationary Landau-Lifshits equation 

MaAM = 0, M^^a{r), (2.1) 

where M = M(a;, y) = {Mi, M2, M^), r = {x,y) e and a{r) is the square of the 

nominal magnetization, which is an arbitrary function of r, in general. 
The equation 

aAM+{Ml + Ml)M^^MAa (2.2) 
is easily obtained from (2.1). We set 

i=l 

where CTj, i = 1, 2, 3, are the Pauli matrices. In terms of the matrix M, Eq. (2.2) is 

+ (aM,M-^), = a,,-/, (2.3) 

where z = x -^iy and / is the 2x2 unit matrix. 

Equation (2.2) (or the equivalent equation (2.3)) describes an inhomogeneous magnet, 
a generalization of the two-dimensional version of the stationary Landau-Lifshits magnet. 
Since the question of the integrability of these equations remains unanswered in the general 
case of an arbitrary function Q;(r)), we consider the simplest version, where Eq. (2.3) is 
integrable, obtained by setting 

4Aq; = = 0. (2.4) 

Then (2.3) implies that a two-dimensional isotropic magnet whose squared saturation 
moment is a harmonic function of the coordinates is described by a nonlinear cr-model in 
a curved space. 

We note an interesting analogy with gravity theory. Equation (2.3) under condition 
(2.4) has the same form as the gravity equations with two commuting Killing vectors for 
the part of the metric g containing the off-diagonal terms [8-11]. The difference is that 
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we have M = M* (where * denotes the Hermitian conjugation) and detM = —a here 
and g = and detg = —o? in the gravity case. 
We set 

3 
i=l 

Therefore, = 1, i.e. S is the normahzed magnetization vector. In accordance with 
(2.3) and (2.4), we now have 

{aS^Sh + (aS.S), = 0, (2.5) 

which is the equation of a nonhnear 0(3) a-model in a curved space. If a = 1, Eq. 
(2.5) becomes the equation of the standard 0(3) cr-model (the two-dimensional isotropic 
Heisenberg ferromagnet), for which the boundary- value problem was solved in [3]- [4]. We 
also note that the above-mentioned case of gravity reduces to a similar 0(2, 1) cr-model 
in a curved space and that dimensional reduction of the axially symmetrical self-duality 
equations in the Yang-Mills theory also leads to a similar model. The analogue of S is 
then given by g*g, where g G SL{N, R) [12]. 

We first consider some simple automodel solutions of (2.5). We take the ansatz 



cosx sinxexp(— i$) 
sinxexp(i$) — cosx 



(2.6) 



where x — xi'^- P)j ^ = ^{oi^P) is a real- value function, and f3 is the harmonic function 
that is conjugate to a in the sense of the Cauchy-Riemann conditions. Substituting (2.6) 
in (2.5) produces the system of equations 

{^aa + (y~^^a + ^^p) Sm X + 2($«Xa + ^/3X^) COS X = 0, 

(2.7) 

2{xaa + a-\a + Xpfs) = (*a + sin 2x 

For X = 0, we have S = a^. For x = 7r/2, the function $ satisfies the Laplace equation; 
it then follows that 

$(q;,/3) = j dR{s)Zo{isa)e''^, 

where Zq{x) is the Bessel function (Jq or Nq) and dR{s) is a spectral measure. Finally, 
ior ^ — kp with a free parameter k, the function x = x(ct) satisfies the equation 

-1 . . 

Xaa + a Xa = y sm 2x, 

which reduces to the third Painlevc equation. 

System (2.7) is analogous to the equations that arise in gravity theory [13]; how- 
ever, one of those equations resembles the Liouville equation rather than the sine- Gordon 
equation as is the case here. 

Equation (2.5) is the compatibility condition for the linear 2x2 matrix system 

aSzS ^ ,^ aSzS ^ 

(dz-^^)^^O, a, + ^^^ = 0, 2.8 
g+ a Q — a 
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where 

Q = g{a, (3,u) = 1(3 - u + {u - {u + 7), 7(z) = a + (2.9) 

7 is the function that is analytic within the sample, and -u G C is a "hidden" spectral 
parameter that does not depend on the coordinate. Thus, Equation (2.5) can be solved by 
a general version of ISTM. This version, however, has not yet been fully described in view 
of the technical complications (some progress was recently made in [14] for the gravity 
theory equations). Therefore, we use the dressing method for deriving exact solutions. 
Because the strategy in this case is essentially the same as in [15], we give only parts of 
the calculations that allow us to obtain the final answer. 

To system (2.8), we add the reduction conditions and the formula for reconstructing 
the potential (magnetization): 

^i-\u) = **(-«), S^{u) = *(t(m)) J, S = *(oo+)C, 

(2.10) 

jr^j^^i, cc*^i, j = c*(oo-), 

where r is the transposition of the sheets of the Riemann surface F of the square root 
entering (2.9) and the signs + and - pertain to the upper and the lower sheet, with the 
upper sheet being fixed by the condition that g ^ eis u ^ 00+. Further, the 2x2 
matrices C and J are introduced to select different versions of the background solutions 
(the classical vacuum). In addition, the involution conditions 

g{u) — —g{—u), g{T{u) — a^/ g{u) (2-11) 

hold. Relations (2.9)-(2.11) allow us to write the answer. We seek the solution of (2.8) in 
the form 

•^{u)^x{u)^%u), (2.12) 

where '^^{u) is the bare solution of (2.8) corresponding to the vacuum and the matrix 
x{u) is to be determined. We further assume that the matrices x(w) and x~^(w) are 
meromorphic functions on the Riemann surface F and have a finite number of poles 
(those of the x matrix occurring at the points u = Vi), 

where /Xj = g{vi), Vi are some complex numbers specified below, and the standard Dirac 
notation for vectors is used. 

Obviously, if x has a pole aX u = Vi, then it has a pole on the second sheet of F at 
u — Vi — T{vi). It follows from the condition x(h)x~^ {u) — I and Eq. (2.13) that 

\m >= E > <h\ = - < N,, - < ^^1^^- > 



Hi - Vj 



Equations (2.8), (2.12)-(2.13) allow us to find 
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where < di\ and |cj > are arbitrary bra and kct vectors. The first of the reductions in 
(2.10) gives < Qil = \pi >* and Ui = —jli, and the second one imphes that < q| J =< c|. 
Here, the vector < q] corresponds to the pole contribution coming from u — Vi and < Cj| 
to the contribution oi u — vi — T{vi). 

The relations obtained allow us to write the " A^-soliton" solution of (2.5); for simplicity, 
we limit ourselves to the " one-soliton" solution 



S = S^ + 



< q\S^\q > 
L»(|//|2 + q;2) 



S^\q > ® < q\So 



< q\q > 



Here, n = q{v) 
constant vector, and 



a^D{ix + //) 



ll\q> ® < q\S^ + iJi.S^\q > 



-v)\c >, < q\ 



\q >* 



\q> ® < q\ 

(2.14) 

® <q\ . 

|c >e is an arbitrary 



D 



a2 



< q\q >^ < q\S^\q >" 



(|/x|2 + a 



2^2 



Equation (2.14) allows one to construct the "N + 1-soliton" solution starting with an 
"A^-soliton" back-ground. 

We next analyze the exphcit form of the solution obtained. The simplest choice for 
the background, 5"° = (73, corresponds to C = /, J — (73, — (73. Setting < c\ — (1, c) 
(without losing generality), we have from (2.14) that 



S+ — Si + iS2 
where 



'1- IcI 



+ (l + |c|^) 



1- 



2|cp 



, (2.15) 



D 



1^ 


2 


[(1 + 


c 


2)2 a\ 


1 - 


(' 








2 + q;2)2 



When c = — iexp(2i$), Eqs. (2.15) are simplified considerably and become 

S+^- sin 2^6^^*, ^3 = - cos 29, 



(2.16) 



where 6 = arg/x. This reproduces ansatz (2.6) with $ = const and x = 7r-|-2^. Moreover, 
it can be shown that x satisfies the second equation in (2.7) with a vanishing right-hand 
side. For a purely imaginary v, it follows from (2.11) and (2.16) that the excitation over 
the -S"" = 1 vacuum is absent. This is true also in the general case described by Eq. (2.14). 
We also note that solutions (2.15) and (2.16) have no nontrivial analogues in the case of 
a constant saturation moment. 

Another possible choice of the bare solution, = exp(— i/c/Sas) ai, where k is a real 
constant, leads to the equations C = J = ai. In this case, system (2.8) can be explicitly 
integrated as 

*°(a,/3,^^) = e'='^^(^('^'^'")-'^^)/^ 

where g is defined in (2.9). Again setting < c| = (1, c), we can generalize the formulas de- 
scribing magneticovortex excitations [16] to the inhomogeneous case under consideration: 
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S+ = e^*^^ + 



4\c\ 



aD 



cos{ri + ia) cos rj cosh ^ 



cosh — i9) 



(2.17) 



where 



2|cp 



D 



4|^|2 



sinh ^ cos 77 — i^^ — ^ cosh ^ sin 77 



cosh ,^ 



(|//|2 + q;2)2_ 

In these formulas, ^ = A;//r + In |c|, 77 = k{iii — (5) — arg c, 9 — arg /i, a — ln{\ii\/a). The 
formulas for 5*+ and 6*3 give the magnetization distribution characterized by local variation 
scale ~ |A; y /i^l and the wave number ~ V if^i ~ Z^)!; which correspond to 
the spatial modulation of the magnetic structure of the sample. For KL ^ 1, relations 
(2.7) describe an inhomogeneous spiral magnet that can produce multiple harmonics of 
the fundamental frequencies. 

The above "solitons" allow us to estimate the energy functional evaluated on these 
solutions. To do so, we assume that the function a depends on the temperature in addition 
to the spatial variables: a = a{r,T). The exact general form of this dependence can be 
found by considering the system on a lattice [17,18], which is beyond the scope of this 
paper. We therefore merely estimate the energy and the heat capacity near the critical 
point. 

According to the Landau theory [19], the functions a and /? behave as follows near 
the Curie point T^. 

a(r,T)^/o(r)t, /3{r,T) ^ go{r)t, (2.18) 

where t = {T^ — T)/Tc. The second estimate in (2.18) follows from the first one and fact 
that the quantities a and P (as well as /o and go) are related by the Cauchy-Riemann 
conditions. Taking a and P to be small and t — > 0+, we see from (2.9) that 



(2.19). 



The system energy is defined as 

E=^ [ dxdy{\/M)'^ = I [ ^^'^^I" 



+ a(vS)' 



where Q is the domain occupied by the sample. This relation can be rewritten as 



E = [ d\[aTr(S,S, 
Jn 



+ 



2a 



El + Eq. 



(2.20) 



(2.21) 



Using 



we obtain 



d^r — dxdy — "^"^7 ^ d'j 



da /\d[3 = ——d'j A d'y, 



1 r dadp f dadp ,^ 

Eq^ - , El = / — — Tr{S^S^), 

8 Jn' a Jn' Iz 



(2.22) 
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where Q' is the domain in the plane of the variables obtained via the conformal 

mapping 7(2;) of fl from the {x,y) plane. It follows from (2.22) that 

^0 ^1{P2- Pi) In — ~ (Tc - T) ln(T, - T), (2.23) 

from which, in particular, we can find the contribution to the heat capacity as 

dE 

Co ^ — ^ ln[e{T, - T)]. (2.24) 

Estimating Ei and Ci = dEi/dT requires knowing the explicit form of the solution for 
S. We first consider the simplest solution (2.16). We then find 



dadp 
7^ 7 

with 9^ and 9^ satisfying the relations 



Ei^8 [ a9,9^ = 8 / dadpa9^9^, (2.25) 

Ja' 7. 75 in' 



^7 = ^(M|))7, % = ^(M|))^ (2.26) 

Using the deformation method developed in [20] , we can show that 

/ 4q; \ 1 1 / 4q; \ 1 1 

+ [t-T^] = —r + T—T^ (2-27) 



ji — aj^ jji — a /X + a y/x + o;/- jji — a fx + a 



from which, setting /i — g, we obtain 



_ 1 g + g . _ 1 g + g 



2i {g + a) {g — a)^ 2i (^ — a) (^ + 0;) 

Thus, in accordance with (2.25), we have 

E^ = 2 [ dad(3 ^^^4^- (2.28) 

Using (2.18) and (2.19) and taking the size of to be finite, we finally obtain the 
following relations satisfied near the Curie point: 

Eo^tl^t, CoRilni, Ei^t^, ci^t^. (2.29) 

Thus, the background heat capacity satisfies the law characteristic of the two-dimensional 
Ising model [17, 18] with the critical index equal to zero, and the energy and the heat 
capacity of the excited state over the vacuum satisfy powerlike laws. 

For convenience in calculating the energy of more complicated solutions, we rewrite 
Eq. (2.21) as 

to/ 1 

E = / dadp [ (2^3.^3.- + l^+.l' + l^+.-n + (2.30) 

J^' 7;^ 7^ 80; 

We consider solution (2.15), setting |//| = atan^ and $ e [0, 7r/2). Then, 
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where 



D = 



c 2|c 
[-Bcos2^ + iAtan9], 83^! 



(2.31) 



4 cos^ 6 



sin^ D = D, 



9 = arg/i, A = 1 + |cp, B = 2 — A. Using the deformation method, we obtain 



cos^ $ 



a 



Q — Q — 2a 



1 



[Q-\-a){Q — a) 2 



tan$ 



(2.32) 



cos^$ 



-7^ 



1^1 



Q — Q — 2a 1 ^ 

7 -r^ r — - tan ^> 

\Q — a){Q-\-a) 2 



Taking (2.31) and (2.32) into account, we obtain 

£^0, El fti t Int, Co, Ci 



Int, 



(2.33) 



which means that the background energy and the heat capacity of (2.15) give the same 
contributions as the corresponding quantities evahiated for the states with excitations. 

We consider solution (2.17). Using the same notation as in the previous example, we 
write it as 



2\c\ 



e**^^ [cos (77 + ia) sin 2$ cos 77 — 



a'^D 

2|cp 
aW 



cosh,^ 
cos 6* 



cosh(e- i9)i 



— cos 2$ sinh ^ cos rj + tan 9 cosh ^ sin 77] . 



Evaluating the corresponding derivatives and estimating the integrals as before, we again 
obtain 

Eo, El ^tl-at; Cq, Ci ^ hit. (2.34) 

The approach described here contains, therefore, three main ingredients: the integra- 
bihty of the model, the dressing method, and the deformation method. Because of its 
universality, the approach can be used with other interesting applications of (integrable) 
two-dimensional models. 

To conclude this section, we show how the transition to the homogeneous case can be 
accomplished in the framework of the inhomogeneous magnet model. We consider system 
(2.8) and expression (2.9), which contains the "hidden" spectral parameter u. Setting 
X — g/a and letting a tend to unity, we obtain 



g = il3 -u + \J{u- 1)(m + l- il3) = A. 



(2.35) 



Because 7(z, z) = a{z, z) + il3{z, z) with z = x + iy is axi analytic function, the Cauchy- 
Riemann conditions imply that ^ /3o as a — >■ 1, where j3o is a constant which we set 
equal to zero. Then we see from (2.35) that u = — 1/2( A + 1/A), and in the language 
of spectral parameters, therefore, the correspondence between the homogeneous and the 
inhomogeneous case is given by the "mirrored" Zhukovskii function. 
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3. The Pohlmeyer transformation and gauge equivalence of 
two-dimensional integrable boundary-value problems 



Wc consider the nonlinear 0(3)(T-model on the half-plane R"^ — {{x,y) : — oo < x < 
+00,1/ >0}: 

S,,- + (S,S,-)S = 0, (3.1) 

where S{x,y) = {Si, S2, S3), |Sp = 1, and = {l/2){Sx — iSy). Model (3.1) coincides 
with the elliptic version of the n-field model [21] and is one of the versions of the chiral 
field model. We assume that S(x, 0) and Sy(a;, 0) are given smooth functions such that 
(unless otherwise stated) 

S(a;, 0) ^ (cos 2a;, sin2a;, 0), ^ 00, (3.2) 

i.e., S belongs to the class of functions describing spiral structures [22]. Following [23], 
we also assume that the relations 

= = 1, (3.3) 
hold on i?^. These relations then imply the constraints 

S^-S^ = 4, S.S, = 0, (3.4) 

which are obviously compatible with (3.1). 

Setting f{z, z) — S^Sz, we perform the Pohlmeyer transformation assuming that 

Szz = ciS + C2Sz + caSg, 

where q = Ci{z,z),i = 1,3, are some complex- valued functions. It can be easily verified 
that ci = —1, C2 = — //z/(l — /^), and C3 = fz/i^ — P)- This implies the equation 

= 1 - - (3.5) 

This equation is integrable, because it is the compatibility condition of the linear matrix 
system 

-^.^A-^, = (3.6) 

where \t' = '^{z,z), A = A{z,z), B = B{z,z) e Mat{2,C), and the matrices A and B 
have the forms 

Here, A e C is the spectral parameter, and |/| > 1. 

In what follows, we use the notion of the "phase" space of a nonlinear elliptic equation, 
which can be introduced by a natural analogy with the hyperbolic case [21]. We do this 
in the example described by Eqs. (3.1) and (3.5). 

The generating functional of the model given by (3.1) and (3.2) is 

= dx dy [-{Si + Si) + -j{S' - 1) + u^SSy - 2J , (3.8) 
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where I'jyj = 1,2, are the Lagrange multiphers. The phase space Ai^^^ consists of the 
smooth functions Si{x, 0) and Siy{x, 0), i = 1, 3, and the analogues the Hamihon equations 
are 



with the " Hamiltonian" of the model given by 



dx 



Si) + 2 



(3.9) 



(3.10) 



The Poisson structures of the phase space is defined by the nonvanishing fundamental 
brackets with the constraints involved in (3.8) duly taken into account [21]: 

{Sayix), Sbix')} = [5ab - Sa{x)Sbix)]5{x - x'), 
{Say{x), Sby{x')} = -[Say{x)Sb{x) - Sby{x)Sa{x)]6{x - X'). 

For arbitrary smooth functionals F and G, we then have 

■ 6F SG 6F 6G 



+ 



' 5G ^ 
5S„ 



[F.G] 



dx< 
5S„ 



6Sy SS 6S 6{6yS) _ 



+ 



+ 



5S, 



'6G^ 



We now consider Eq. (3.5) assuming that f{x, 0) and fy{x, 0) are given smooth real- 
valued functions such that f{x,0) — > 1 and fx{x,0), fy{x,0) — > 0, as |x| oo, with the 
second term in (3.5) vanishing at infinity. The generating functional of this equation is 



dx 



1 f+f 

- ^" ^^ +4(1-/) 



2 P-1 

The "phase" space A^*^-^) is spanned by the functions q = ln|/ + y/p — 1| and p 
fyl \I\P — ll, I/I > 1, and the equations of "motion" are 

qy = {H^f\q}, qyy = {H^f\qy}, 
where the "Hamiltonian" is given by 



dx 



1 if!) 

2/2- 



4(1-/) 



and the Poisson structure of the "phase" space is defined by the fundamental brackets 
{p(x), q{x')} = 6{x — x').^ 

Equation (3.5) is useful because it generates a broad family of integrable elliptic equa- 
tions. In what follows, we assume (unless otherwise stated) that u = u{z, z) is a smooth 
real-valued function defined on and, in addition, that 



u{x, y) -^0, |a;| 



0. 



(3.11) 



■^We stress that it is not our aim here to construct the Hamihonian formahsm for eUiptic equations. 
In this case, unfortunately, we do not have "intuitively obvious" variables of the angle-action type, and 
all the analogies with the hyperbolic variables become rather conventional. 
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Setting f{z, z) — coshw, we obtain the equation 

A'u = -4sinhM (3.12) 

from (3.5). Equation (3.12) emerges in the theory of two-dimensional Boltzmann-Poisson 
systems at negative temperatures [24] and also in some reductions of the Chern-Simons 
model [25]. It has a geometric meaning, describing the embedding of a negative-curvature 
surface in three-dimensional Euclidean space [26] . 

Using the GE of the problems given by (3.1), (3.2) and by (3.11), (3.12), which was 
proved in [4], [5], and also using Eq.(3.4) and the definition of /, we obtain 

S^ = 4cosh^|, Sj = 4sinh^|. (3.13) 

These equations, which relate the "potentials" of the two boundary- value problems, imply, 
in particular, that the GE does not impose any constraints on the corresponding "phase" 
spaces (provided, of course, that / > 1). 

We list some other possible choices of the functions f{z^ z) and the corresponding 
equations that lead to interesting, important physical applications. 







— cosh u, 


ddu — sinh u, 


(3.14) 




z) = e±". 


ddu — - 


dudu 
-2 sinh 14 ± — — , 

g±2« _ I ' 


(3.15) 


f{z, 


z) = -e±" 


, ddu — 


dudu 
2 sinh u ± , 

g±2« _ I 


(3.16) 




f{z,z)^ 


^ cos u, 


ddu — ± sin u, 


(3.17) 




f{z,z) = 


± sin u, 


ddu = ± cos u. 


(3.18) 



We consider the question of the GE of the problems given by (3.1) and (3.11), (3.17). 
For this purpose, we note that Eq. (3.1) can be rewritten as 

CikiSkddSi = 0, 

where cm = e^, with eiki being the totally antisymmetric tensor in R^. Therefore, Eq. 
(3.12), as well as Eqs. (3.14)-(3.16) corresponds to the model of a compact-manifold 
magnet. This situation changes for Eqs. (3.16)-(3.18), and the sought for GE is achieved 
on a noncompact manifold of matrices S e SU (1,1). We set 

The condition det^ = —1 and Eq. (3.19) imply that 

Geometrically, this means that the end of the vector S = {S^, S"^, S^) moves over the 
surface of the one-sheet hyperboloid (in contrast to the unit sphere in the "compact" case). 
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(3.19) 



The matrix S can be expanded with respect to the basis elements of the Lie algebra su{l,l) 
as S = S'lii, where tTc^tt^ = rjap + iclpT^^, = 1,2,3, r] = {r^^^} = diag{l, 1,-1) 

is the KiUing metric, ^ are the structure constants of su{l, 1), and the commutation 
relations are [7ro,7r/3] = 2ic^^ vr^. 

If Eq. (3.17) is now taken with the + sign, it follows from (3.4) and from the expression 
for / that 

{Sir + {Sir-{S!f = 4sin'^, 

(Sir + {Sir- {Sir = -^cos' I 

and the requirement for conservation of the asymptotic behavior implies that the GE 
to (3.1), (3.2) is achieved under the condition that u{x,0) — > 7r{mod2'7i). Then, the 
constraints on the space Ai^^^ of the corresponding 0(1, 2) a-model take the form of the 
system of differential inequalities 

< {Sir + {Sir - {Sir < 4, 
-4 < {Sir + {Sir - {Sir < 0. 

The GE conditions for the other equations in the above list are derived in essentially 
the same way. 

Equations (3.15) and (3.16) seem new. It would be interesting to find applications of 
these equations, which could be considered the Boltzmann-Poisson equations (for positive 
temperatures in the case of (3.16) and negative ones in the case of (3.15)) generalized to 
the case of a quantum-classical system, and also to quantize them. We do not rule out 
other possible applications of these equations in physics and geometry (a certain regular 
method to obtain (2 + 0)-dimensional equations in the differential geometry of surface is 
given in [27])[|. 

The elliptic version of the relativistic field theory model proposed in [28] takes the 
following form in the one-component case: 

^p- , 2n wdwdw _ 

oow = —w(l — w ) —■, w = w (3.20) 

1 ~ 

Equation (3.20) can be naturally called the elliptic Getmanov model (a similar equation 
for a complex- valued function was considered in [29]). It is known that in the hyperbolic 
version, this model is gauge equivalent to the sine-Gordon equation in Minkowski space 
[28]. The solution of (3.20) can be easily related to the solution of Eq. (3.5) as / = l — 2w'^. 
This allows us, in particular, to obtain from F^-^^ the density of the generating functional 
for (3.20) as 

F(^) = 8{dw dw/{w^ - 1) + w^). 

Properly speaking, the instances of GE established above between a number of mod- 
els on the half-plane and the model described by (3.1) and (3.2) require a more rigorous 
(mathematical) proof, which should include the demonstration of transitions of these 
models to the cx-model. Clearly, this demonstration can be done in each particular case 

■^Exact solutions of Eq. (3.16) were constructed in [7] using the Darboux transformation method; the 
relation of these equations to the Bitsadze equation was also considered there. 
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using the appropriate associated linear problem and constructing the gauge transforma- 
tion matrix. 

4. Exact solutions of the sinh-Gordon equation 

We consider a plasma consisting of electrons and singly charged ions. We assume that 
the electron and ion temperatures are the same and each plasma component has relaxed 
to the Boltzmann distribution. In this case, the dimensionless potential of the electric 
field u — e$/T satisfies the Boltzmann- Poisson equation [30] 

Aw = 4sinhw, (4.1) 

where A is the two-dimensional Laplace operator written in dimensionless variables, r = 
R/2rD, e is the electron charge, T is the temperature, $ is the potential, R is the radius 
vector, ro — T / {SnenQY^'^ is the Debye radius, and Uq — = noi are the electron and 
ion concentrations. 

Equation (4.1) also emerges in other physical problems: in the theory of strong elec- 
trolytes and conducting films [31], in the 0(2, 1) cr-models [32], in the calculation of the 
electrostatic contribution to the DNA molecule free energy [33], and so on. The same 
equation is obtained by the dimensional reduction (from four to two dimensions) of the 
self-duality equations for the S'f/(2)-valued Yang-Mills fields under a special choice of 
the ansatz [34]. It was assumed in [34] that a spontaneous gauge-symmetry breaking 
occurs: part of the potential components tend to constant matrices as r — > oo with half 
of the them playing the role of Higgs bosons thereby leading to the massive field theory 
corresponding to (4.1). 

We also note that there arc sufficiently rigorous procedures to derive (4.1) from the 
Bogoliubov-Born-Green-Kirkwood-Ivon chain of equations [35]. 

We assume that (4.1) is defined on the half-plane i?^, and that 

u(x,y)^0, r = ^ oo, (4.2) 

where u — u{x, y) is a sufficiently smooth real- valued function. 

It is useful to consider first the linear version of (4.1) and (4.2) that corresponds to the 
Debye-Hiickel approximation. In this case, the solution obtained in the Fourier-transform 
parametrization is 

u{x,y) = - r'^hsiX) e^^W^-'W^ 

TT JO A 

(4.3) 

hB{\) = i y dxe-'^^^^'' [l{\)u{x, 0) - Uy{x, 0)] , 

where k{X) = \ — and /(A) = A + A^^. It is assumed that the surface charge density 
at the boundary Uy{x, 0) decreases faster than exp{—2\x\) as — >• oo. Then 6s(A) can be 
analytically continued from the real axis to the domain in the plane of the complex variable 
X — Xr + iXi that is bounded by the strophoid curve A|j(2 — A/) — A/(A7 — 1)^ = and 
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that includes the unit semicircle in the right half-plane. In this domain, 6b(A) = 6b(1/A) 
because u is real. 

Solution (4.3) occurs under the boundary conditions 



u 



6b(A) = 0, A = A<0, 

[x, 0) + - / dx'Ko{2\x - x'\)uy{x', 0) = 0, 

TT J 



(4.4) 
(4.5) 



which follow from the requirement that -u be bounded asr ^ oo. In (4.5), Ko{x) is the 
Macdonald function. Equation (4.5) is the Fourier transform of (4.4). The asymptotic 
form of (4.3) as r — > cxd found using the saddle-point method is 



u{r,a) ~ (nr) ■2bB[iexp{—ia)]exp{—2r), 



(4.6) 



where x = r cos a and y = r sin a. Expression (4.6) describes the effect of the linear Debye 
screening. 

We return to the original problem formulated in (4.1) and (4.2) and assume that the 
given functions u{x, 0) and Uy{x, 0) cannot be arbitrary in a given class of functions (for 
example, in the Schwartz class) and are constrained by some condition (analogous to 
(4.5)) to be formulated in what follows. 

The auxiliary linear system of equations corresponding to (4.1) is 



. A cosh u , 



Uz 

}(^3 - ycr2 - 



sinhw 
'2\ 



-ai 



,A coshw .Uz sinhii 

-^2 + ^)^^-^T^^-^^^ 



(4.7) 



(4.8) 



We introduce the matrix of Jost solutions (4.7), which are determined by the asymp- 
totic formulas 

^± = (^± ^J) = ^ ^^^^^^ ^ ^ ^4 9^ 

We set 



, A = A, 



(4.10) 



where T{\) is the transition matrix. It has the properties following from (4.7) and (4.10) 



r(-A) = <72T(A)<72, T(A)=T(^), (4.11) 
and also the unimodularity property detT(A) — 1. From (4.11), we have 

a(A)=a(i), 6(A) =6(i), (4.12) 

where the second involution takes place, strictly speaking, for A = A. It, however, admits 
an analytic continuation into some domain in the complex plane (we do not need more 
detailed information about that domain in what follows). We set (jr^ = _ 
(0f , 02^). Then the Volterra integral equations for the direct problem, which are equivalent 
to (4.7), become 
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<pf{x, A) =ej + j dx'gf{x - x' , A) Q{x' , X)<pf{x', A), (4.13) 
where j — 1, 2, ci = (1, 0)^, 62 = (0, 1)^, the Q matrix is 

( cosh u—l sinh u i iu^ \ 
2a 2A 2 1 
sinhM iuz, coshn— 1 I ) 
2A 2 2iA / 

and gf are the bare Green's functions 

gi{x, A) = Td{Tx) diag{l, e'"'''), 

(4.14) 

g2{x, A) = ±^(±x) diagie''''', 1) 

It follows from (4.13) that 

(t)t, 02 G i/(i^+), 0r, 0^ e i/(i^"), (4.15) 

where H{Q) is the class of functions that are analytic in the domain fl, and Q"*" = {A : 
ImX > 0}. 

Taking (4.10) and (4.13) into account, we have the integral representations for the 
scattering data a (A) 6(A), 

a(A) ^1- j dx<el, Q{x, A)(/)^ (x. A) >, (4.16) 

6(A) = / dxe-'^'^^^'' < ef, g(a;, A)0+(a:, A) >, (4.17) 



where <, > denotes the scalar product of vectors in C^. In the standard way, Eq. (4.8) 
implies the relations that describe the spectral data evolution, 

a(|/,A)=a(|/,0), 6(|/,A) = 6(0,A)e-'W^ (4.18) 

from which it can be easily seen that the condition that 6(|/, A) be finite for any y > 
implies 

6(A) = 0, A = A < 0, (4.19) 

and a gap thus appears in the spectrum of the associated linear problem (see also [3, 
36]. Turning to (4.17), we see that a condition similar to (4.5) emerges in terms of the 
functions n(,T,0) and Mj,(,t, 0). 

Wc consider in more detail the properties of the coefficient a(A). It follows from (4.16) 
that a(A) = 1 + 0(1/A) as |A| 00 and a(A) e H{Q+) and the unimodularity of T(A) 
and Eq. (4.19) imply that a(A)a(— A) = 1. Therefore, using (4.16), we conclude that 
a(0) = 1. Assuming now that there exist simple zeros of the coefficient a(A),a(A„) = 
0, ImXn > 0, n = 1, ...N, and recalling the condition u — u together with Eq. (4.12), we 
obtain one more representation for a(A), 

"W^n^ n ' 7 ' Yy (4.20) 
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where 2Ni + 2N2 = N and the zeros A„, 1 < n < 2Ni, belong to the unit circle (the 
analogues of kinks for the hyperbolic version of the sin-Gordon equation [21]). The zeros 
Xjn and 1/Xjn: 2iVi + 1 < m < 2 A'"! + N2 constitute an inversion with respect to the unit 
circle (the analogues of breathers [21]). It is obvious that there is no topological charge 
in the model, in contrast to the equation Au = sin-u (see [37, 38]). 

It follows from the first relation in (4.18) that Ina(A), in particular, can be considered 
a generating functional of the integrals of "motion". The standard procedure brings the 
system of equations for the elements of the column (f)^ to the Riccati equation 

+ QuF^ = -ik(X)F + (Q22 - Qn)F + Q21, (4.21) 
where F — 02i/0ii- Further, we obtain 

Ina(A) = - rdx[Qu + Q12F]. (4.22) 

Setting 



F = ^F_fcA^ |AHO, F = ^FfcA-^ |A| 

k=0 k=l 



00, 



00, 



R — Qll + Q12F — ^2 ^-k^'' ^ I'^l 0, — ^ i^fcA |A| 

k=0 k=l 

we see from (4.21) that 

R^i = i[{coshu-l)/2 + ul/2]+iul/2 - idl\n{coshu/2) - {coshu - 1) /2d^{u^/ sinhu), 

i?o = «9xlncoshM/2, Ri = -i[{coshu~l)/2 + ul/2), R2 = -{l/8)d^iul), 

and so on. Since the left-hand side of (4.22) does not depend on y, the rigth-hand side 
can be evaluated at y = 0, which gives an infinite number of conservation laws. 
Now let 

00 

Ina(A) = ^IsX^\ \X\ 00 

s=l 

and 

00 

Ina(A) = ^ /_,A^ |A| ^0. 

s=0 

It follows from (4.12) that /, = /_„ hs = and Iq = 0. Using (4.20)-(4.22), we arrive at 
the trace identity 



^[EAr^+ E {X'r' + K''^')]- fdxR2s-rix), (4.23) 

m=l n=2Ni+l •' 

where s = 1, 2, .... The functions i?2s are either equal to zero or given by total derivatives 
in X, their integrals vanishing in view of (4.2). 
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Wc now solve the inverse problem. The standard argument employing the analytic 
properties of solutions leads to the Riemann problem of reconstructing a piecewise- 
analytic vector function from its jump on the boundary (A = A): 

^i^;^ = 02 A) + r{y,X)e''^'^^<j>t{x,y,X), (4.24) 

where r{y,X) = b{y,X)/a{X) = 6(0, A)/a(A)e-^(^)^ = r(0, A)e-'(^)^ is the reflection coeffi- 
cient. After some calculations, we finally obtain the system of singular integral inverse 
scattering equations 

^^^^^^^^ I rj^^^ 

~, A + A„ Jo zn u + X + lU 



n=l 



«„ = e:+.i;!!^^<^..A(^,y)- r| (4.26) 

„=1 Xm + Xn Jo Zn IJ, + Xm 

where = 2(A^i-|-iV2) and m„ = hn/a\Xn) are the discrete spectrum transition coefficients 
such that rUn — —X^fhn for 1 < n < 2A^i and m„ = — A^m^+jva for 2A^i + 1 < n < 
2A^i + N2. In Eqs. (4.25) and (4.26), further, 4>i{x, y, A), and 0i'„(x, y) are the respective 
eigenfunctions of the continuous and the discrete spectra. 

The formulas for reconstructing the potential should be added to the system given by 
(4.25) and (4.26). Setting 

00 

0+(x,y,A) = ^0+(x,y)A^ \X\^0, 

k=0 

we insert this expansion into the matrix equation as 0+ and compare it with the cor- 
responding expansion (4.7). We thus obtain one of the formulas for reconstructing the 
potential. 



zsmh- = 2. m„(y) - / — (l)Ji{x,y, n)e . 4.27) 

2 Xn Jo 2m n 

A similar expansion for |A| — > 00 yields 

^^^^ = E ^t.ne''^'"^^rnM ' f ^r{y, /.)0+ (x, y, ,.)e^^^^^\ (4.28) 

Relations (4.25)-(4.28) allow us to construct the simplest solutions of Eq. (4.1). We 
first consider the "no-soliton" case a(A) = 1. Then system (4.25), in which r — > 00, can 
be easily iterated, the first iteration giving ^^(r. A) = 1 -|- 0{l/r). Using (4.27) and also 
the fact that u{r, a) — > as r — > 00 for a e [0, tt], we obtain 

Jo TT /I 

where /(/x) = ik{n) cos a — /(/x) sin a. Assuming that 6(/i) can be analytically continued 
into the required domain in the complex plane, we obtain 
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(4.30) 



Comparing this with (4.6), we see that the only difference in the two expressions is the 
replacement &b — > 6 and that (4.30) becomes (4.6) with sufficiently weak fields. This 
effect can be called a quasi-linear Debye screening: the scale of the screening is the same 
as in the linear case, and the influence of strong fields amounts to renormalization of the 
preexponential factor - the effective charge of the boundary - which contributes to the 
asymptotic behavior of the potential. 

We now set 6(A) = 0, < A < oo, which means that we consider the case of reflection- 
less potentials. The system of inverse scattering transform equations then reduces to a 
system of linear algebraic equations. Recalling (4.28), we obtain the " A'"-soliton" solution 



u{x,y)^2\n—, A± = det(5^„ ± exp ^-[fc(A„) - ^AJ]). (4.31) 

where 1 < m,n < N . One of the simplest solutions (the analogue of a breather) corre- 
sponds to setting Ni = 0, N2 = 1, Ai = I/A2 = exp(7 -|- iO), and mi(0) — A^m2(0) = 
exp((5 -|- ip). It follows from (4.31) that 



, , ^ + h , sin[2r sinh7cos(^^ + a) + p — 6*)! coth7 / , „„x 

u{x,y)^2hi -. h = ■ ^ — . (4.32) 

^ ^' l-h^ cosh[2rcosh7sin(^ + a) +7-(5-ln^2i^] ^ ^ 

In the general case where {9 + a ^ 0, tt) , this solution describes a charge-density wave, 
which falls off as r — > 00 according to the law 

M ~ exp(— 2r cosh7| sin(6' + q;)|. (4.33) 

The parameters involved in (4.32) can be related to the values of the fleld at the boundary 
y = by the trace identities (4.23), which now become 

cosh7COS^ = — — f dxux{x,0)uy{x,0), (4.34) 

coshjsine ^ J dx[4:{coshu{x,0) -1)+ ^{ul{x,0) - ul{x,0))]. (4.35) 

As can be seen from (4.32), the field does not fall off for + o; = tt, which means that 
a long-range order occurs in that direction, i.e., there is a coherent structure describing 
the distribution of the electric field and plasma component densities with a considerable 
charge separation. We also note that for h — ±1, solution (4.32) becomes singular (we do 
not consider this case here). 

Another simple solution (the analogue of a kink) that follows from (4.31) is obtained 
by taking A^i ^1, 0, A„ = e^^", n = 1, 2, mi,2(0) = i e^^-^+'^^'^ : 

cosh[2r sin(^) cos(^ + a) - ^] cot ^ 
sinh[2r cos(^) sin(^ + a) - ^ - ln(i tan ^)] " ^ ^ 
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This case corresponds to the nonUnear Debye screening. As r — > oo , we have 

h ^ ^-2r.Sr.[e,^a)+25, ^ ^-2. sin(e,+a)+ 25, q ^ ^1+^ + a < TT. (4.37) 

The term with the smallest coefficient in front of 2r survives in the exponent of (4.37). 
Because this coefficient is always less than one, the field falls off more slowly under the 
nonlinear Debye screening than in the quasi-linear case, and, moreover, the fall-off is 
anisotropic. The parameters characterizing the fall-off (4.37) can be related to the bound- 
ary values of the fields using trace formulas (4.34) and (4.35) if we replace cosh 7 there 
by cos ^1^, and also Q by 

The slower fall-off of the potential under the nonlinear Debye screening can be ex- 
plained by the pressure of the electric field, which "pushes aside" the electrons and ions 
in the polarization cloud surrounding the boundary. The collective nonlinear phenomena 
observed here - the threadlike "radiographic examination" of the plasma with the electric 
field and its slower fall-off law - seem to take place only in the two-dimensional case. The 
field between the piontlike charges then decreases more slowly than in three dimensions, 
which explains the appearance of long-range correlations. 

We now consider several thermodynamic relations and also those involving energy. 
The full energy of the medium (using dimensionless quantities) is 

E = y da;rf7/[(VM)^ + 4MsinhM]. (4.38) 
In the linear case, this relation can be represent as 

E = ^ dxdy[{Vuf + 4u^] = Ef + Ef + Ep. (4.39) 
Using (4.3) and the relation 

J (i^e^[fe(A)+MM)]a. ^ (27r/(A l{X)))S{ii - 1/A), 

(x) (v) 

we evaluate Ej^ , Ey and Ep with the result 



where 



1{X) 
2 /-"ociA 



E = -l ^|6b(A)|2. (4.41) 

TT Jo A 

Restricting ourselves to the "no-soliton" sector in the nonlinear case, we can exactly 

(x) (v) 

evaluate the quantities analogous to those found above only for the sum Ey + Ey and 
only for the continuous spectrum (the integrals diverge in the other cases). From (4.28), 
we see that 

Ux tUy dfx , \ 4. / \ 



4 JO 27riyU 

where we introduce the function xf (x, y, A) = 0f(a;,|/, A)e'^W^\ It now follows that 

/ dxdy{ul -ruD^A ^r(0, /.)r-(0. A) dye-^^^ ^^^^^^ J dxxUx, fi)xUx, A). 

(4.42) 
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To evaluate the integral, therefore, we must know the orthogonality relations for the 
functions Xii{^: The method for obtaining these is to write the integral equations for 
Xii{x, A) and Xii(^) A) and then multiply them using the formulas 

J l{fj,) k{X) — k{fj,) 

where P denotes the principal value and Tq = e-''''^^'^ ii+ik{n) > 5; e^'=(^)«2-^'=(^)«i, 

S > ; 1) — ^2- After some calculations, we find 

and finally 

/«i,(«^ + «;) = ^f + (g^ + l). (4.43) 

In the linear limit as m — > (a(/x) = 1) and — > 6s (/^), this expression becomes equal 
to (4.41). 

We consider the problem of the charge of the medium (as a whole). Its total value is 
4 / dxdy sinhu ^ / dxdyAu ^ - / d^Uy(^,0). (4.44) 



Whenever the solution u is singular, a term given by the sum of integrals along the contour 
that cuts out the singularities should be added to the right-hand side of (4.44). Using the 
"breather" solution and evaluating this sum, we can show that this term vanishes and we 
have (4.44) again. This relation, therefore, shows that the electric neutrality condition 
should be understood such that the total charge of the medium and of the boundary 
vanishes. 

The trace identities (4.34) and (4.35) also have a very clear physical interpretation. 
To show this, we introduce the Maxwell stress tensor of the electric field: Tq,^ = E^E-y — 
(l/2)5a^E^, where E = {Ex,Ey) is the electric field vector. Because the field has a 
potential, we obtain 



where 



[curlE , E]« = d^Ta-y - 2To a«(cosh — - 1), (4.45) 
Tj, = $5$„ T,, = -h^-<^l), (4.46) 



2 

^ = {VT/a)x, rj = {VT/a)y and a = V^ttc^^. It follows from (4.45) that 



Tr,r,-2noT{cosh^-l) =-d^Tr,^. (4.47) 
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On the other hand, evaluating the force apphed to the boundary from the medium, we 
obtain = 0, and 

F„ = 2eno [ d^dr]smh = -27^0 / ^^(cosh ^ - 1). (4.48) 

Jr'^ T J T 

Comparing (4.35)-(4.48) with (4.34) and (4.35), we see that the "trace identities" involve 
the quantities that can be expressed through the Maxwell stress tensor and the force 
applied to the boundary. 

We consider several properties of the parameter 7, entering (4.34) and (4.35), where we 
assume that 7 = 7(T, no). We note that the "breather" solution degenerates as 7 — > 0, all 
the eigenvalues of the auxiliary linear problem becoming equal. This degeneration occur 
in two ways: () at uq — ric, where T is arbitrary; (b) at T = Tc where uq is arbitrary. 
Expanding the right-hand sides of these identities in the power series in (n—ric) or (T—Tc), 
we find the behavior of 7 near the critical parameters ric and Tc. 

In case a, it follows from (4.34) that 7 ~ \j2{nc — n)/nc, as n /" ric- This behavior of 
7 occurs under the additional condition, which follows from (4.35), 



/ rfe(cosh^ - 1) = 0, = $(^^^e,0). (4.49) 

A similar analysis in case b results in 7 ~ \J'2{Tc — T)/T etsT y Tc, with condition (4.49) 
replaced by 

4 / d^(cosh ^ - 1) = / sinh ^, = $(-^, 0). (4.50) 

The degeneracy of the eigenvalues can probably be interpreted in this case as a certain 
"phase" transition in the nonlinear medium. This hypothesis, however, requires a further 
analysis. 

Note also, that the equation 

Au = -4sinhii, (4-51) 
describing the plasmas at the negative temperatures, was solved in [39] . 
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